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Abstract. Using a phase field crystal model we study the structure and 
dynamics of a drop of colloidal suspension during evaporation of the solvent. 
We model an experimental system where contact line pinning of the drop on the 
substrate is non-existent. Under such carefully controlled conditions, evaporation 
of the drop produces an ordered or disordered arrangement of the colloidal residue 
depending on the initial average density of solute and the drying rate. We obtain 
a non-equilibrium phase boundary showing amorphous and crystalline phases 
of single component and binary mixtures of colloidal particles in the density- 
drying rate plane. While single component colloids order in the two dimensional 
triangular lattice, a symmetric binary mixture of mutually repulsive particles 
orders in a three sub-lattice order where two of the sub-lattices of the triangular 
lattice are occupied by the two species of particles with the third sub-lattice 
vacant. 



1. Introduction 

The phase field crystal (PFC) technique has been used extensively to study a variety of 
phenomena like solidification, dendrite formation in binary alloys [1], order-disorder 
transitions, and the elastic [2j |3] and rheological response of solids under external 
loads. The technique off'ers the unique advantage of being able to model atomistic 
detail within a coarse-grained, free energy functional approach. While there have been 
attempts to relate PFC [H |3] models to more realistic dynamical density functional 
theory [5J [7] approaches, the basic simplicity of PFC is appealing. 

In this paper we present our results on the application of the PFC model in its 
simplest version to a complex problem of great technological interest viz. assembly of 
nano-particles from a drying suspension |51f^ll0j. The experiment consists of placing a 
drop of suspension containing colloidal particles on a substrate and allowing the solvent 
to evaporate . Under suitable conditions, which we explain shortly, the residue may 
consist of particles arranged in a regular periodic order. Such regular arrangements of 



nano-particles have many applications from functional arrays of sensors to templates 
for the production of more complicated nano-structures. 

The physics of the drying process is quite complicated and may be briefly 
described as follows. The rate of evaporation at the contact line, where the drop meets 
the substrate, is much higher than that of at the middle of the drop due to higher 
curvature [12j . So, the evaporation of the drop proceeds mainly from the contact line 
[13j . As the solvent dries, a radially outward flow of solvent is set up which replenishes 
the evaporated liquid [21 [TSl [TS] and tends to advect particles from the center to 
the periphery. Simultaneously the contact line moves inward towards the center of 
the drop with evaporation, dragging particles along with it due to surface tension 
forces. These competing flows result in a buildup of solute at the contact line which 
effectively pins it. Subsequently, with the stationary contact line, fluid evaporation 
changes contact angle [TJ such that the pinning force decreases suddenly and the 
contact line becomes unpinned starting the cycle again. Such repeated pinning-de 
pinning transitions leave concentric rings of solute in the wake of the drying drop 
leading to the so-called "coffee ring" effect. Inhomogeneities of the solute density 
produced by this process leads to a disordered residue with the fluid having to flow 
through a highly ramified network of stationary and moving particles. Analysis of this 
complete process is complex and to our knowledge a complete solution to this problem 
does not exist. Also the disordered residue patterns produced makes it impossible to 
use this process as an effective technique of nano-particle assembly. 

From the above discussion, it is clear that suppressing contact line pinning would 
result in a drastic reduction of complexity and has the promise of producing an ordered 
residue [TH [TB] . Indeed, recent experiments on the drying of colloidal suspensions on 
carefully prepared substrates do show such drying induced ordering [TTJ (THl [20] . A 
mobile contact line which moves uniformly leads to a steadily shrinking drop where 
solute particles experience only an inward compressing force arising from the moving 
meniscus. The radial compaction leads to solidification once the density surpasses 
the freezing density. However, if the drying rate is too large, the liquid evaporates 
before this density is reached and the colloidal particles are stranded in a disordered 
arrangement. 

We have, in the past, modeled this system of uniformly drying colloidal suspension 
on a substrate without contact line pinning using Brownian dynamics simulations 
of particles interacting with a coarse-grained solvent field [TT]. Here we use a PFC 
approach to simulate the system using two coupled coarse-grained fields corresponding 
to the solute and the solvent local densities respectively. Preliminary results of this 
study has been published as a conference proceedings [3T]. Finally, we have also 
generalized our system to represent a symmetric binary mixture of particles and we 
consider the case where the particles repel each other. In this case, there are three 
coupled coarse-grained fields representing the local densities of the two species of 
particles and the solvent. 

The paper is organized as follows. In the next section, we introduce our PFC 
model for single component suspension. This is followed by details of the numerical 
method used to solve the resulting coupled dynamical equations. In section 4, we 
present our results for the single component system. We generalize to the binary 
system in section 5 and present the results. Finally we end the paper with a brief 
discussion of our results and pointing out future directions of further work. 
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2. Phase field crystal model for a single component suspension 

We begin by considering the evaporation of the solvent as a two dimensional (2d) 
liquid-gas phase transition. The phase separating fluid is described by a free energy 
functional !Fi containing quadratic and quartic terms and an external chemical 
potential as follows: 

= j [-a'/'' + + c(V0)2 + hdpl (1) 

Here (/)(r, t) is a scalar order parameter field (local density of the solvent) at position 
r and time t with a, b and c > as phenomcnological parameters. The free energy 
functional here is the standard '0*' model of 'liquid-gas' phase transition with the 
external field '/i', representing a chemical potential which favors < (vapour) regions 
i.e., drives the evaporation process. The time-evolution of the fluid is determined by 
the time dependent Ginzburg-Landau (model A) dynamics [32], 



= r(2a0-4503-H2c(V^(/))-/i) (2) 

The drying rate of the solvent can be controlled by varying the parameters h and 
r, the kinetic coefficient. For small values of h the free energy density in Eqjl] has 
minima at (pi and 02 where (pi > Q and (/)2 < 0. We denote (j) — (pi as liquid phase 
and = 02 as vapor phase of the solvent. For /i > 0, the solvent starts to evaporate 
and when all the solvent dries up then </> becomes 02 everywhere. 

For solute particles we consider a PFC model in its simplest form defined by a 
dimensionless free energy functional [21 [3] is 

F,m]- Jd'r{t[e+il + vr]i' + ^} (3) 

where is a coarse-grained field whose extrema corresponds to the equilibrium 
positions of the atoms and therefore may be interpreted as a dimensionless density 
field [2] of the solute. The phase diagram of this system, obtained by minimizing Eq[3j 
has been extensively studied in Ref. [2] . A schematic of the phase diagram obtained in 
Ref . [2] is reproduced in Figjl] The stable phase for all e > is the uniform phase where 
tp = constant = ipav, where ipav is the spatial average of ip. As e becomes negative, 
the free energy is minimized by modulated phases with modulations in either one or 
two dimensions. Further, the phase diagram in the ipav — e plane is symmetric in ipav 

There are some inherent difficulties involved in using the PFC model to describe 
the drying induced ordering of colloids. Firstly, there is no evidence for a smectic 
phase in a colloidal system with spherical particles, so this phase is unphysical and 
needs to be discarded. Secondly, unlike the usual freezing phase diagram where a solid 
is obtained from the liquid on cooling or increasing the density, in the PFC model, 
the transition is re-entrant so that there are limits to how much ipav can be increased 
at fixed e such that one obtains a periodic solid from the uniform phase. Finally, 
the order parameter jump across the freezing transition can become very small for 
small |e|, this is in variance with a real liquid - solid transition which is strongly first 
order and remains over all physically accessible temperatures and pressures for most 
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Figure 1. A schematic of the phase diagram in the Van — e plane from Ref. [2] 
obtained by minimizing the the Landau free energy functional T2 (Eq. |3} over 
choices of functions V'Ci") ■ Note that the phase diagram is symmetric along ^av = 
reflecting the basic symmetry of the functional. There are three distinct phases 
where (1) 1/) is uniform as in a liquid, (2) is periodic only in one direction (smectic 
or stripe phase) and (3) has the periodicity of the two dimensional triangular 
lattice. The shaded areas are coexistence regions and the arrows show quench 
paths connecting liquid and triangular solid phases. 



common substances. In spite of these caveats, it is possible to work within a range of 
parameters such that meaningful results are obtained as we describe below. 

In our case as the solvent starts to evaporate continually from the contact line, the 
liquid-gas interface shrinks and consequently the solute particles are forced to come 
together along the radially inward direction and hence, the value of ipav increases. 
The evaporation of the solute therefore induces a force on the particles coupling the 
fields (f) and tp in the dimensionless free energy J-"2[{^/'}]. Rewriting ^2[{V'}] including 
an appropriate (leading order) coupling gives us 

-^2 = |d^r{|[.+ (l + V^m+^+0^} (4) 

where the coupling coefficient has been chosen to be unity. 

The equation of motion for conserved field tp (using model B dynamics pflHl]) is 
therefore given by, 

^ = ^.(W^) = + (1 + V^)^]V. + + 0) } (5) 

In the above equation fi is the mobility of the solute which is governed by the 
following equation. 

^i = (tanh(7<?!>) + l)/2 (6) 

The functional form of fi should encode the fact that if solute particles leave 
the solvent during evaporation, their motion ceases due to the lack of Brownian 
fluctuations in gas where = (/)2 < 0. So we choose the functional form of fj, such 
that as 4> = (f>2, /i — >■ 0. The value of 7 controls the rate at which /i vanishes in the 
limit (j) (f>2, near the liquid- vapour interface. Here the average density of the solute 
particles and the drying rate are the control parameters. The average density of the 
solute particles can be varied by varying ipav The drying rate can be controlled by 
varying h and F both. 
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As the hquid- vapour interface shrinks continuously enhancing ipav, it eventually 
reaches a range where a periodic ip is expected and an ordered (triangular) phase forms 
(see Fig{T]). On the other hand, if the drying rate and the initial value of is such 
that the final value of tpav remains below the critical value needed for solidification till 
the final stage of evaporation, then a disordered solute phase is expected. 

In this model two important time scales are responsible for the 'drying induced 
ordering' phenomena. The first, ti, is related to the drying time of the solvent and 
other, T2, is related to the relaxation of the solute particles to its crystalline minima. 
No ordering in solute particle is possible if tl^av is lower than the critical value and 
Ti < T2, because the solute particle does not get the sufficient time to reach the 
crystalline minima in that case. Ordering is only possible if ti > T2 and the ipav is 
greater than the critical value. 

We now give below details of the procedure used to solve the dynamics for the 
model given in Eqsj2] and [5] The coupled dynamics of the solute and the solvent are 
not exactly solvable by analytical methods in general. However, the dynamics of a 
drying drop of pure solvent without solute is solvable using a reasonable ansatz for 
the shape of the drop as we show below. We utilize this fact to define two methods of 
solution. In method I [3T], we solve for the dynamics of the drop radius analytically, 
which is then used as an input to solve the dynamics of the particles, i.e. i/"- In 
method II, the dynamics for both the coupled fields (j) and are solved numerically. 
Both methods gives similar results but method I is faster than method II. Below we 
will describe both in order. 

In method I |21| . for a given ^av, we use an ansatz for (f) remembering that it 
represents a two dimensional drop of solvent; 

(/.(r,i = 0) =01 for\r\^R 

^0^+ (^l_i^(i?_r) for R<\r\<{R + 

= (j)2 otherwise (7) 
Substituting this profile in the free energy Eq[2| and minimizing it with respect 
to the interfacial width, we obtain ^ = 2^ c(4>i-ih y^ for -| <C 1, with Xs ~ 
I --a(0|-0|) ^ 6(0?_0|) ^ 1^ ^j^g ^^^g ^j^^^gg ^^j.^^^^ ^^^j^ obtained 

as, 



d4> dR d(j) 
dt dt dr 



r=R 

which, after some straightforward algebra, gives the dynamical equation for the 
radius of the drop during evaporation as 

+ M*.i-,f2)] + | (9) 

Where we have assumed that evaporation depends only on radial coordinate such 
that the shape of the drop is kept fixed and ^ 1. After integrating Eqjoj we get 



where the Rq is the initial radius of the drop and D — (0^-^3)2 ^ ^2) + ^(^^i ~ 

</.|) + ;i(0i-02)]. 

Since Eq{lO] is a transcendental equation, we solve the Eq{lO] numerically using 
the bisection method at every time step. Figj2] shows typical solutions of R{t). In this 
simulation we choose ip{v,t — 0) randomly [23] over space with a non-zero average 
value. At every time step we get the value of the radius of the drop solving Eq[TO| 
numerically and with this value of radius a new configuration of solvent density 0(r, t) 
is generated according to Eq[7]at every time step. Using this value of 0(r, t), we update 
tp{r,t) using Eq. [5] A Fast Fourier Transformation along with Euler discretisation is 
used to solve Eqjs] When the value of the radius of the drop becomes equal to R* 
such that ^ < R* << Rq, we get the final configuration of ijj and the time required to 
arrive at this stage, starting from Rq is the total drying time (i^j). 




Figure 2. (a) and (b) show the variation of the radius of the solvent (R) with 
time for = 1568.15, h = 0.006 (fast drying) and = 6620.32, ?t = 0.002 (slow 
drying) respectively. R* = 50.0 in the both cases. Fig. is taken from 21 . 



Method II, on the other hand, is more straight-forward, where we solve for 
both solvent as well as solute dynamics numerically, again using an explicit Euler 
discretization scheme with a suitably small time step. At each time step we update the 
density profile of the solvent and use the updated profile to solve the solute-dynamics 
as before. 

Finally, we fix the length and time scales for simulation and choose a set of values 
for the parameters of the model as follows. The length scale is set by the coefficient 
of the squared gradient term in the free energy for ip and the time scale is set by the 
coefficient of the time derivative of ip, both of which are set to unity. We take 7 = 10.0 
in Eqj6] to set the mobility of the solute particles inside the liquid as well as outside 
the liquid. The parameters a = 0.008, b — 0.004, c = 1.0. The value of h is adjusted to 
control the drying time. Using these values, i.e., values of a, b, c, we get the minima of 
the dimensionless free energy J"i at (pi = 0.651, (j)2 = —1.15 for h = 0.006 (a typical 
value of h that is used here for simulating fast drying) and (pi = 0.930, 02 = —1.06 
for h = 0.002( a typical value of h that is used here for simulating slow drying). 
The interfacial width becomes, ^ = 31.77783, for h = 0.006 and ^ = 44.65193, for 
h — 0.002. A system of size 512Ax x 512Aa; with Aa; = 1.0 is taken with periodic 
boundary conditions. Time step dt = 0.01. We have taken the radius of the drop 
as Rq = 200 when time t = and R* = 50 when time t = td- Here F = 1.0 and 
e = -0.375. 
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3. Results and discussions 



The final configuration of the solute depends on its final value of ipav ■ As mentioned 
before, if it is within the range of the values of the phase diagram [5] where we get 
a periodic triangular solid, the resulting structure will be ordered. If it is below the 
range we get liquid-like disordered state and otherwise we get the stripe phase, which 
is unphysical for our context. We typically begin the drying process with an initial 
low solute density or small value of ipav, and this value and the drying rate together 
determines the final structure. 




(a) X (b) X (c) 



Figure 3. (a),(b) and (c) are the final configuration of i/'(r,t) for ip^^ = —0.7, 
= —1.0 and ■i/'J^j, = —1.2 respectively. All the above three figures [21) are for 
td = 1568.15, /i = 0.006 and F = 1.0 and they are showing three phases, striped 
phase, ordered (triangular) and disordered liquid-like phases respectively. (Using 
Method I) 




Figure 4. (a) final configuration of i/'(r,t) for ipav = — 1.2,td = 1568.15, /i = 
0.006 and F = 1.0 (fast drying) and (b) final configuration of ip{r,t) for tpav = 
-1.2, id = 6620.32, ft = 0.002 and F = 1.0 (slow drying). For same tpav but for 
different drying rate (a) shows liquid-like disordered phase and (b) shows ordered 
(tringular) phase |21| . (Using Method I) 



We show typical final configurations of ip for various initial density and drying 
rate. Both the methods described earlier, produce similar results. While Fig 3 a) 
shows a striped pattern in the final configuration of the solute, Figjsjjb) and Fig 3 'c) 
show ordered and liquid-like final configurations respectively. Here the initial values 
of tpav are taken according to the order: ip^^ > tjj^^ > tp^^ where the superscripts 
denote respective figures. The total drying time (t^) in all these cases are same. 
This implies that for a fixed drying rate and for small initial densities, they are in 
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liquid-like disordered phase and for large initial densities they are in stripe phase. 
We get crystalline phase in the intermediate initial densities. Figjl] (a) shows a 
typical disordered liquid-like configuration for fast drying and Figj4] (b) shows a typical 
ordered (triangular phase) configuration for slow drying where both have tpav — —1.2 
initialy. So, for same initial solute density we get crystalline order when the drying 
process is slow and liquid-like disordered structure when the drying process is fast. 
The above results are obtained using method I for solving the dynamical equations. 

When method II is used, the results are quite similar validating our ansatz 
for cj). In Figjsjwe show typical results showing ordered (Figjsf^a)) and disordered 
structures (FiglSib)) produced by varying ipav We have also shown the structure 
factors corresponding to these configurations in Figjsjc) and (d) respectively. We 
have calculated the structure factors ('5'(k)) using, 

S{k) = { / i/;(r)cos(k.r) dr} + { / sin(k.r) dr} (11) 





Figure 5. (a) and (b) are the final configuration o{ip(r,t) for ip^^ = —1.1 and 
= —1.3 respectively. AH the above figures are for same td,h,T which are 
td = 1505.8, /i = 0.006 and F = 1.0. For same drying rate but for different tpav 
(a) shows the ordered phase and (b) shows the disordered phase, (c) and (d) are 
the structure factor of (a) and (b) respectively. (Using Method II) 



In Fig[6] similar results showing ordered and disordered final configurations (and 
the corresponding structure factors) obtained by varying the drying rate keeping ipav 
fixed are shown. Note that both F and h has been varied here to change drying time. 
We have seen that the results depend only on the drying rate and not on h or T 
independently. 

From the previous discussions, we understand that drying rate of the solvent and 
initial density of the solute (both of which are externally controllable parameters of 
the problem) infiuences the final structure of the solute through a dynamical order- 
disorder transition. Here we concentrate on the order-disorder phase boundary of the 
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Figure 6. (a) final configuration of ip(r,t) for ipav = — 1-3, t^j = 1505.8, /i = 
0.006 and F = 1.0 (fast drying) and (b) final configuration of ij){r,t) for ijiav = 
-1.3, td = 1810.26, h = 0.002 and V = 3.0 (slow drying). For difl'erent drying rate 
but for same il}av liquid-like disordered configuration is found in (a) and ordered 
(triangular) configuration is found in (b). (c) and (d) are the structure factor of 
(a) and (b) respectively. (Using Method II) 




Figure 7. Phase boundary of order-disorder transition for two different h i.e., 
drying rate, in the parameter space of total drying time of the solvent and the 
average density of the solute. 



problem. As we are interested only in the two dimensional order-disorder transition, 
we do not allow the dynamics to take the system into the stripe phase, which is 
unphysical to the present context. To ensure this, we take recourse to the following 
artifice. We make /.t not only dependent on 0, but also on ■0 mimicking the slowing 
down of dynamics in the solid phase. We set a cutoff value for ip above which /i = 0.0 
and below which /i satisfies equation([6|. This ensures that the value of ij^av is small 
enough so that at the end of drying, we have only crystalline or liquid-like disordered 
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state, and not the stripe phase. When the total mass of the solvent i.e., integrated 
value of (j) over space, is less than a preset cut-off value representing the detectability 
limit for solvent, we stop the drying simulation and record the final configuration of 
V'(r, i). The time required to arrive at the final stage from the initial configuration, 
is the total drying time {td)- Averaging over many independent initial configurations, 
finally we obtain the order-disorder phase boundary in the parameter space of initial 
solute density and total drying time, which is shown in Figj?] To construct the phase 
boundary, we have used several drying rates obtained by considering two different 
values of chemical potential: h — 0.002,0.006 and four different values of F and for 
each of these values, we have also considered several initial ipav values of the solute. 
The shape of the phase boundary, shown in Fig[7j shows that as we increase the solute 
density, less time is required for drying in order to get the final ordered state. 

4. Generalization to a binary mixture 

In this section we generalize the PFC model of drying induced ordering to a symmetric 
binary mixture of two species of colloidal particles which interact with one another. 
We aim to see under what conditions an ordered crystalline arrangement of the colloids 
is obtained. Needless to say, an ordered arrangement of two species of colloids has 
many possible technological applications [55]. 

A binary mixture can be modeled by a free energy with minimal coupling of the 
form [3] 

F,^ jd\[^[e + {l + V2)^]V.i + ^ + + (1 + ^'f]^2 

+ + (^(V-i + V'2)} (12) 

Here the parameter a couples the two fields ipi and ip2 representing the two species 
of colloidal particles. Note that in Eq{T2]the two fields are treated on the same footing 
and the free energy is unchanged under a switch between the two f/'iS. For a > 0, 
the system prefers 4'i and il)2 to peak at different locations which leads to a decrease 
of the overall free energy. Positive a therefore favors ordering. On the other hand 
if a is negative, the two species of particles tends to occupy, statistically, the same 
positions within the lattice and therefore denotes a disordered crystal where the two 
species of particles randomly occupy lattice sites with equal probability - a randomly 
substituted binary alloy. 

In order to determine the ground states of the above free energy for the (infinite 
size) binary system we consider the following ansatz for ipi and V'2- 

(i) (A) Randomly substituted triangular lattice: This is the structure closest to 
the single component case where the two species of particles occupy sites of a 
triangular lattice randomly with equal probability. In this case ipi and ip2 have 
identical forms given by: 

= V2 = At[cos{qtx) cos(gty/\/3) - cos(2gty/%/3)/2] + Vat, (13) 

(ii) (B) Two sublattice ordered square lattice: This structure consists of a two 
dimensional square lattice which is bipartite and hence can be decomposed into 
two sublattices. One sublattice is occupied by Vi whereas the other is occupied 
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by 11)2- The corresponding forms for ipi and 1/^2 are: 
= As[cos[qsx) + cos(gsy)] + ^av 

= As[cos{q,{x + a/2)) + cos{q,{y + a/2))] + ^^v (14) 

(iii) ( C) Triangular lattice with three sublattice order. It is difficult to make an ordered 
arrangement of two species of particles on a triangular lattice which is tripartite 
due to geometric frustration. In this case, the problem is solved by keeping one 
of the sublattices vacant and ipi occupying the other two leading to a relatively 
open honeycomb like structure. The functional forms of ipi and ip2 are given by: 

tpi = Ast[cos{qstx) cos{qsty/V3) - cos{2qsty/V3)/2] + ipav 

"02 = Ast[cos{qst{x + a/2)) cos{qst {y + a/2V3)/V^) 

- cos{2qst{y + a/2V3)/V3)/2] + ^av (15) 

In the above three equations At, As, A^t are unknown constants and qt — qs = 
qst — 27r/a where a is the lattice constant. 




a a 

Figure 8. (a) Phase diagram in the ipav — ol plane of the three different 
lattice structures. The red line with filled triangle is the phase boundary between 
randomly substituted triangular lattice and triangular lattice with one sublattice 
vacant and the blue line with filled square is that between triangular lattice with 
one sublattice vacant and two sublattice ordered square lattice, (b) plot of free 
energies of the three different structures per unit area as a function of a for 
i/)ai> = —0.25. Green, red and blue colour lines represent randomly substituted 
triangular lattice, triangular lattice with one sublattice vacant and two sublattice 
ordered square lattice respectively. 



Substituting Eq 13 into Eq |12| and after integrating over unit cell and minimising 
with respect to At and qt, we get free energy per unit area of the randomly substituted 
triangular lattice as 

At^ 



45 3 

\ I I \ ^ \ \ 2 I rat) 

+ (e + a + l)T/)au H — 



where At — 



-36V'a 



15(e + a) 



„ _ \/3 

9i - ^ 



'3(e + a)] 

(16) 

and S is the area of the unit 



cell. Similarly, substituting Eq{T4] into Eq[T2j free energy per unit area of the two 
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Figure 9. (a) Plot of the maxima of and ip2 showing the three sub- 
lattice order. Colors yellow and blue corresponds to the maxima of 4'1 and -02 
respectively. Here average value of ipi and ip2 = —0.6, F = 1.0, h = 0.006 and 
total drying time = 1450.0 (b) Comparison between phase boundaries of order- 
disorder transition of single component system and binary system for h = 0.006 in 
the parameter space of total drying time (td) of the solvent and average density of 
the solute. For binary system, i/)^^ is the intial average density of one component 
of the binary system (i.e. ipi or ^2)- 



sublattice ordered square lattice can be written as 



S 



= -A/ + As^[3^Jav^ + e~a] + {e + a + Ijt/jav^ + 



(17) 



where As — \\/ —i'fpav^ — e + a, = 1 and finally, substituting Eq 15 into Eq, 
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free energy per unit area of the triangular lattice with one sublattice vacant becomes 

Eli 
s 



^Ast' - l^^avAst' + ^[9^a.' + 3(e - a/2)] 



+ ie + a + 1)V^.' + ^ 



where Ast — 



-36^Pav^ - 15(e - a/2) 



- V3 
2 • 



In Fig 



(18) 

b) we have 



plotted the free energies of the three structures per unit area as a function of the 
parameter a for = —0.25. It is reflected from the FigjsJ^b) that for a < 0, 
randomly substituted triangular lattice is more probable than the other structures 
because of its lower value of free energy but for < a < 0.3 triangular lattice with one 
sublattice vacant has lower free energy value and for a > 0.3 two sublattice ordered 
square lattice has the smaller free energy value. Figjsjja) shows the phase diagram of 
the three different lattice structures in ^av — ct plane. It is obtained by calculating free 
energies of different lattices per unit area for different values of "0011 and a i.e., putting 
different values of tpav and a in EqjTBl EqfT7| and Eq{T8] A global minimization of 
the free energy among the three chosen structures then yields the phase diagram. 
Figjsja) also informs that randomly substituted triangular lattice is more probable 
in the region where a < and for slightly higher positive value of a two sublattice 
ordered square lattice has the greater chance to appear and in between these ranges 
of a, triangular lattice with one sublattice vacant is more probable. 
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Having determined the phase diagram for the binary system, we are now in a 
position to address the problem of drying induced ordering of the binary mixture. 
The dynamical equations for ■0^ may be written down in conformity with the single 
component case as, 



In our case, a = 0.1 is taken. We also considered here a binary system having 
equal average density of each component. And we found that when each component 
formed an ordered configuration then the density peak of each component does not 
coincide each other in space which is clearly seen in the Figjoja). For this value of a, 
the final configuration is that of an ordered triangular lattice with one sublattice vacant 
(C). We have also tried with higher value of a (positive), where a two sublattice ordered 
square lattice is expected. However, in this regime we always obtain a disordered 
configuration and no equilibrium square lattice. We suspect that this is because the 
circular boundary in our system is more compatible with a triangular rather than a 
square lattice, a situation compounded by the (necessarily) small system size. 

The dynamical phase boundary of the order-disorder transition of binary system 
in the parameter space of total drying time of the solvent and total average density 
of the solute is shown in Fig|9][b) which also indicates, like single component system, 
that as we increase the average density of solute, less time is required for drying to get 
the ordered state finally. We know that total drying time can be varied by changing 
the value of chemical potential h and mobility F but here total drying time of solvent 
is varied by only changing mobility F, keeping chemical potential h fixed {h = 0.006). 
We have also Plotted the phase boundary of the order-disorder transition of a single 
component system in the same figure (Figjojjb)) showing that, qualitatively, the two 
systems are similar to each other. 

5. Discussion and future directions 

In the real svstem[ll). the dynamical phase boundary is not monotonic as in our 
studies (Figsj?] and |9]^b)), but re-entrant signifying that for higher densities one 
obtains, again, a disordered final state. The reason for this, as explained in Ref.fTT] 
is that there are two timescales in the problem. The first timescale corresponds to 
the time required to compactify the system and this decreases with increasing density. 
The second timescale corresponds to the nucleation time for the system to find the low 
free energy ordered state in configuration space. The second timescale increases with 
density due to an impending glass transition. In the simple versions of the PFC model 
that we have used, there is no glassy solution present and hence we cannot model this 
second phenomenon within our approach in any simple manner. There has however 
been some work in trying to generalize the PFC model to include glassy states [Ml US] 
and slow dynamics. We hope to extend our work in this direction in the near future. 

Apart from spherical particles (whether monodisperse or bidisperse) one may try 
to study the drying of non-spherical particles by suitably modifying the PFC model. 
In this case one may need to invoke additional order parameters such as the local 
director. Similarly, in order to incorporate hydrodynamic effects, additional fields 




(19) 




(20) 
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corresponding to the momentum densities of the solute and the solvent need to be 
included. Work along these lines are in progress. 

Acknowledgements. Funding from the Indo-EU project MONAMI is gratefully 
acknowledged. 
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